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ABSTRACT 

Context, For galaxies hosting supermassive black holes (SMBHs), it has been observed that the mass of the central black hole (M.) 
tightly correlates with the effective or central velocity dispersion (<r) of the host galaxy. The origin of this M. - a scaling relation 
is assumed to lie in the merging history of the galaxies but many open questions about its origin and the behavior in different mass 
ranges still need to be addressed. 

Aims. The goal of this work is to study the black-hole scaling relations for low black-hole masses, where the regime of intermediate- 
mass black holes (IMBHs) in globular clusters (GCs) is entered. 

Methods. We collect all existing reports of dynamical black-hole measurements in globular clusters, providing black-hole masses 
or upper limits for 14 candidates. We plot the black-hole masses versus different cluster parameters including total mass, velocity 
dispersion, concentration and half-mass radius. We search for trends and test the correlations in order to quantify their significance 
using a set of different statistical approaches. For correlations showing a large significance we perform a linear fit, accounting for 
uncertainties and upper limits. 

Results. We find a clear correlation between the mass of the IMBH and the velocity dispersion of the globular cluster. As expected, 
the total mass of the globular cluster then also correlates with the mass of the IMBH. While the slope of the M.-cr correlation differs 
strongly from the one observed for SMBHs, the other scaling relations M. - M tot , and M.-L are similar to the correlations in galaxies. 
Significant correlations of black-hole mass with other cluster properties were not found in the present sample. 

Key words, black hole physics - stars: kinematics and dynamics 
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1. Introduction one parameter power law M. - <x and M. - L relations, 

there is now some evidence for large scatter or upward cur- 
O ; 1.1 ■ The low-mass end of the M.-cr relation vature trend for SMB H masses at the hig h est mass regime 

Ti't, •• i r , u u * i u, i (Ma ~ 10 10 M^ e.g. |G ebhardt etaI]l201li iMcConnell et al.1 
L| , Empirical scaling relations between the masses ol nuclear black om i . I a — r h t""l 701 7 s ) 

. _ ■ holes and properties of their host galaxies in terms of to- 

S><^ tal luminosity L, bulge mass Mbuiee (Kormendv & Richstone 

CI '. Il995t iMarconi & Hunll2003l: Earing & Rix 2004) or the effec- 

. Cy _■ tive (proj ected) velocit y dispersion cr dF errarese & Merritt 2000; 

Gebhardtetalj |20004 |Ferrarese & Ford| |2005| |Gultekin et alj Given this puzzling behavior of the M.-cr relation at high 



P ... 

|2009) disclose a strong connection between the formation his- masses , the question arises if different laws apply for the low- 

tory of supermassive black holes (SMBHs) and galaxies. These mass en d as well. Attem pts to measure lower black-hole masses 

scaling relations can be explained by gas accreting black holes j n dwarf galaxie s (e.g. lFilippenko & Holl2003l:lBarth et al.ll2004t 

which drive powerful jets and outflows into the surrounding in- fXiao et alJl20Tl sparsely populate the relation down to a black- 

terstellar medium, thus affecting the star formation rate and reg- hole mass of ~ 10 5 but not lower. The lower the black-hole 

ulating the further supply of matter onto the SMBH ( Silk & Rees masSi t h e fainter they appear in radio and X-ray emission and 

1 998; |DiMatteoetag|2005fc |Springel et alj2005|; |Hopkins et alj me weaker is the kinematic signature and the smaller the radius 

2007|). Alternative explanations reproduce these empirical cor- f influence. Therefore, detecting black-hole masses lower than 



relations by large numbers of merger events in which extreme iq 5 M is challenging. Reaching for the mass range which nowa- 

values of M./M bulge are averaged out. days is ass i g ned to intermediate-mass black holes (10 2 - 10 5 M , 

While galaxies and central SMBHs in the most common IMBHs) requires extending the search towards different stellar 

mass range, M, - 1O 6 ~ 9 M , seem to closely follow the systems. 
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1.2. Detecting IMBHs in globular clusters 

The possible existence of IMBHs in low-mass 


stellar 


sys- 


more investigations this cluster no longer appears as a strong 
IMBH hosting candidate (e.g. Dull et al. 1997; Baumgardt et al. 
2003, 2005; van den Bosch et al. 2006) and the previously de- 



terns such as globular clusters (GCs) was first suggested by 
ISilk & Aronsl d 19751) while studying X-ray sources in a large 
sample of globular clusters. They concluded that the observed 
X-ray fluxes could only be explained by mass accretion onto a 
100 - 1000M G central black hole. This discovery triggered a 
burst of black-hole hunting in globular clusters using X-ray and 
radio emission but also photometric and kinematic signatures. 

Due to the small amount of gas and dust in globular 
clusters, the accretion efficiency of a potential black hole at 
the center is expected to be low. This makes the detection 
of IMBHs at the centers of globula r clusters through X-ray 
and radio emissions very challenging (Miller & Hamilton 2002; 
Maccarone & Servillat 2008). However, several attempts were 
made to detect radio and X-ray emission of gas in the cen- 
tral re gions and to pro vide a black-hole mass estimate (e.g . 
Maccarone et aT] 120051: lUlvestad et all 120071: iBash et all 120081: 
Cse het all 1201*0 ). Recently. IStrader et al.1 (1201 2|) testedseveral 



Galactic globular clusters for the presence of possible IMBHs 
by investigating radio and X-ray emission. They only found up- 
per limits of the order of 1O 2 M . However, the authors had to 
make various assumptions about the gas accretion process such 
as gas distribution, accretion efficiency, and transformation of 
X-ray fluxes to black-hole masses in order to derive those lim- 
its. Besides X-ray and radio emission, the central kinematics in 
globular clusters can reveal possible IMBHs. 

1.3. The first IMBH candidates in globular clusters 



Driven by the results of ISilk & Aroilsl (ll975l) . iBahcall & Woll 
dl976l) claimed the detection of an IMBH in M15 by measuring 
its light profile and comparing it to dynamical m odels. The first 
claim from kinematic measurements was made by Pet erson et al.l 
( 1989) who found a strong rise in the velocity dispersion profile 
of_M15. Later. iGebhardt et al.l (fl99 7, 2000b) and Ge rssen et ail 
d2002l) estimated a black-hole mass in M15 of (3.2±2.2)xl0 3 Mn 
from photometric and kinematic observations. However, after 



tected X-ray sources turned out to be a large number of low-mass 
X-ray binaries, but new detections of IMBH candidates in other 
clusters followed. 

Using integrated light near the center of the M31 cluster 
Gl. IGebhardt et al.1 d2002[ 120051) measured the velocity disper- 
sion of this cluster and argued for the presence of a (1.8 ± 
0.5) x 1O 4 M IMBH. Furthermore, X-ray and radio emission 
were detected at the cluster center, consistent with a black 
hole of the same ma ss dPoolev & Rappaport 2006; Kong 2007; 
Ulvestad et al. 200 7J). This r e sult, h owever, was recently chal- 
lenged by Mill er- Jones et al.l d2012l) who found no radio signa- 
ture at the center of Gl when repeating the observations. 

Another good candidate for hosting an IMBH at its cen- 
ter is the massive globu lar cluster co Centauri (NGC 5139). 
iNovola et al.l ([2008, 20 1 0) measured the velocity-dispersion pro- 
file with an integral-field unit and used dynamical models to an- 
alyze the data. Due to the distinct rise in the velocity-dispersion 
profile they claim a black-ho le mass of M. — 40 000 Mp. The 
same object was studied by Anderson & van der Mare l (2010) 
using proper motions from HST images. They found less com- 
pelling evidence for a central black hole, but more importantly, 
they found a location for the center that differs from previ- 
ous measurements. Both Gl and co Centauri have been sug- 
gested to be stripped nuclei of dwar f galaxies (Free manlll993L 
iMevlan et al. 2001; Jala li et al.ll2012l) and therefore may not be 
the best representatives of globular clusters. 

Further evidence for the existence of IMBHs is pro- 
vided by the discovery of ultra luminous X-ray (ULX) 
sources at no n-nuclear locations in starburst galaxies (e.g. 
Fabbiano 1989; Colbert & Mushotzkvl [19991: iMatsumoto et al.l 
2001; Fa bbiano et al.ll2.Q0lh . The brightest of these compact ob- 
jects (with L ~ 10 41 ergs -1 ) imply masses larger than 10 3 M o 
assuming no beaming of the X-ray emission and accretion at 
the Eddington limit. Probably one of the best IMBH candidates, 
the ULX source HLX-1, was discover ed in an off-center posi- 
tion of the spiral galaxy ESQ243-49 bv lFarrellet"aT1(l2009l) and 
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iGodet et all (120091) and is most likely associated with a group of 
young stars remaining from the core of a stripped dwarf galaxy 
(e.g. lSoria et al.ll20 1 0ll20 1 3l) . 

In this paper we take the latest results for IMBH measure- 
ments in globular clusters from the literature and search for 
correlations with properties of their host cluster. For the first 
time, enough datapoints on IMBHs are available in order to per- 
form this analysis. Established relations for supermassive black 
holes and their host galaxies, such as the M. - cr relation or the 
M, - M, ot , are the focus of this work. Probing the low-mass end 
of these correlations will provide important information about 
the origin and growth of supermassive black holes. In Section|2] 
we introduce the sample and the data from the literature used in 
this work. Section [3] describes the methods we apply for deter- 
mining correlation coefficients and regression parameters and in 
Section |4] we discuss the major correlations and their statistical 
significance. Finally, Section summarizes our work and lists 
our conclusions. 



2. The sample 

The sample we chose for our analysis consists of 14 Galactic 
globular clusters for which we have kinematic measurements 
of their central regions. The majority of these measurements 
come fro m integral-field spectroscopy and wa s carried out by 
our group ( Liitzgen dorf et al.l20 1 1 [ 120121 120 1 3| Feldmeieret al., 
2013 submitted). 

These clusters were observed with the GIRAFFE spec- 
trograph of the FLAMES (Fiber Large Array Multi Element 
Spectrograph, Pasquini e t al. 2002) instrument at the Very Large 
Telescope (VLT) using the ARGUS mode (Large Integral Field 
Unit). The velocity-dispersion profile was obtained by combin- 
ing the spectra in radial bins centered around the adopted pho- 
tometric center and measuring the broadening of the lines using 
a non parametric line-of-sight-velocity-distribution fitting algo- 
rithm. For larger radii, the kinematic profiles were completed 
with radial velocity data from the literature, when available. In 
addition to the spectroscopic data, HST photometry was used to 
obtain the star catalogs, the photometric center of the cluster and 
its surface brightness profile. For each cluster photometry and 
spectroscopy were combined in order to apply analytic Jeans 
models to the data. The surface brightness profile was used to 
obtain a model velocity-dispersion profile which was fit to the 
data by assuming different black-hole masses and M/Ly pro- 
files. The final black-hole masses were obtained from a^- 2 fit to 
the kinematic data. Table[T]lists all the clusters of the sample and 
their major properties. 

In four clusters of our sample we found signatures of an 
intermediate-mass black hole. These clusters are NGC 6388, 
NGC 6266, NGC 1904 and NGC 5286. All of them show a 
rise in the velocity dispersion and require black-hole masses 
between 1.5 x 10 3 - 2 x 1O 4 M . A high radial anisotropy in 
the center of these clusters would also explain the rise in the 
velocity-dispersion profile. However, N-body simulations have 
shown that any anisotropy in the center of a globular clus- 
ter would be smoothed out after a couple of relaxation times 
dLtitzgendorf et al.ll2.Ql ll) . Since all these clusters are older than 
at least 5-10 relaxation times, we conclude that the best explana- 
tion for the rise of the kinematic profile is indeed a black hole. 
Nevertheless, the result of the cluster NCG 1904 has to be treated 
carefully since the rise might occur from a mismatch of outer and 
inner kinematics. The rest of the clusters in our sample shows 
less evidence for a central black hole with shallow or even drop- 



ping velocity-dispersion profiles. For those we adopt the lcr er- 
rors of the measurement as the upper limit. 

a) Centauri, Gl and NGC 6715 are globular clusters with 
measured black-hole masses not obtaine d by our gro u p. For 
co Centauri we use the latest result of iNovola et al.l ([2010) 
for the black-hole mass of the IMBH in a> Centauri of ~ 
47 000 M Q , which was derived by studying VLT/FLAMES 
integral-field spectroscopy data of the inner kinematics, simi- 
lar to the method described above. This value was also con- 
firmed through a comparison with N-body simulations by 
Ualali et al.1 d2012l) . The globular cluster Gl in M31 was ob- 
serv ed with the Sp a ce Te lescope Imaging Spectrograph (STIS) 
by iGebhardt et al.1 (2002) and the velocity-dispersion profile 
compared with dynamical Schwarzschild models. For our anal- 
ysis we use their derived mass of M. - (1.8 + 0.5) x 10 4 M o . 
Ibata et al.l d2009h measured the inner kinematics of the globular 
cluster NGC 67 1 5 (M54) using a combination of VLT/FLAMES 
MEDUSA and ARGUS data. They found a cusp in the in- 
ner velocity-dispersion profile consistent with a ~ 9 400 M 
intermediate-mass black hole. Unfortunately they did not pro- 
vide any uncertainties on this measurement since they could not 
exclude a high amount of radial anisotropy in the center which 
would als o explain the rise of the profile. With our results ob- 
tained in Liitzgend orf et al.l (1201 ll) however, we are confident 
that high radial anisotropy is not very likely to reside long in 
the centers of globular clusters. We therefore adopt the IMBH 
measurement of llbata et al.l (|2009) and adopt a conservative un- 
certainty of 50% on the black-hole mass. We stress that high 
anisotropy in a globular cluster is unlikely for undisturbed clus- 
ters, however considering the violent environment of some clus- 
ters in the sample (i.e. NGC 6388 and NGC 6266), disk and 
bulge shocking cannot be excluded. The uncertainty in the black- 
hole mass measurement might therefore be underestimated in 
some cases. 

The last two clusters included in our sample are M15 and 
47 Tuc. Both of them seem to have no IMBHs at their cen- 
ters. M15, which was observed with several instruments over 
the last years, was one of the most popular objects for the 
study of IMBHs. As a post-core collapse cluster is a bad can- 
didate for hosting an IMBH at its center since a massive black 
hole would enlarge the core of the cluster and prevent cor e col- 
lapse (iBaumgardt et al.ll2003HNovola & Baumgardtl20Tll) . Also 
47 Tuc is not a good candidate for hosting a massive black hole 
at its center. McLaughlin et al. (2006) obtained HST proper mo- 
tions for 47 Tuc and found no evidence for a rise in its velocity 
dispersion profile, i.e. an IMBH. 

The final sample of globular clusters and their properties are 
listed in Table Q] The cluster parameters were taken from dif- 
ferent sources (see references in Table Q3 and their black-hole 
masses from the observations listed above. 



3. Dealing with censored data 

With our current dataset we encounter several challenges: 1) 
The dataset is small. With n — 14 we enter statistical regimes 
(n < 30) where many statistical approaches, i.e. the Spearman 
correlation test, are not reliable anymore. 2) The dataset contains 
not only measurements, but also a large number of upper limits. 
The statistical method which deals with th ese so called "cen- 
sored" data sets is called survi val analysis. Feigelson & Nelson 
(1985) and llsobe et al.l (1986) proposed several techniques for 
applying survival analysis to astronomical datasets. Especially 
treating bivariate data (a dataset with two variables) with up- 
per limits and determining correlation coefficients and linear re- 
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Fig. 1. M, - cr (left), M. - Mjot (middle) and M. - Ljot (right) relation for globular clusters with IMBHs. The gray arrows indicate 
upper limits. We compare the slopes of the best-fits from four different methods. 

Table 2. Correlation parameters and significance calculated through various methods. 
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gressions is described in detail in these papers. However, sur- 
vival analysis does not treat uncertainties of the detected mea- 
surements. This would bias our result since our dataset 3) also 
contains large asymmetric uncertainties in both variables. 

We are not aware of a heuristic method which treats all of 
these caveats properly. Therefore, we will present several at- 
tempts to analyze the data by accounting for each problem sep- 
arately and discussing their differences. The goal is to test for 
possible correlations between two parameters of the dataset. For 
those where the correlation is significant we obtain the linear 
regression parameters a, (3 and en, where 

y-a+[3x (1) 

and eo indicates the intrinsic scatter i.e., dispersion in y due to 
the objects themselves rather than to measurement errors. 



the correlation coefficient r becomes zero. Positive and negative 
signs of the Kendall-r imply a correlation and anticorrelation of 
the data, respectively. Table|2]lists the Kendall-r and the signifi- 
cance of its deviation from zero (P) for several data combinations 
as calculated from the IDL routine R.CORRELATE. 

As suggested from Table [2] and Figure [TJ there are three cor- 
relations observed in our dataset: M, - cr, M. - M lol and M. - L 
as observed for SMBHs in massive galaxies. In order to compare 
the correlation properties with measurements from galaxies we 
perform a power-law fit to the data defined as: 



log lo (M./M ) = a a +# T log 10 [cr/(200kms- 1 )] 
log lo (M./M ) = a M + J M log lo [M, of /(lO 11 M )] 
log 1() (M./M ) = a L +J3 L log lo [L, of /(lO 11 L )] 



(2) 
(3) 
(4) 



3.1. Partly treating uncertainties 

Since there are no correlation coefficients known which would 
include uncertainties of the measurements, we are limited to 
the standard statistical methods. Furthermore, due to our small 
dataset, we cannot perform a meaningful Spearman correlation 
test. This leaves us with the the Kendall-r rank correlation coef- 
ficient as a measurement for the significance of the correlation. 
The r parameter is calculated through the numbers of concordant 
(ranks of both elements agree) and discordant (ranks of both el- 
ements differ) datapairs and the total number of elements in the 
dataset. If the data x and y are independent from each other, then 



All three relations are fitted using the algorithm described by 
iTremaine et al.l (I2002I) . where the quantity 



X' 



-Z 



[log 10 (M. ;/ ) - a -/3log m (cr/o- )] 2 



i=\ 



et + e 2 



+ fi 



2*2 



(5) 



is minimized and the 68% confidence intervals of a and /J are 
determined by the range for which Ay 2 — x 1 ~Xmm ^ 1- The un- 
certainties of M. (cm) and cr (e^) are considered asymmetrically 
in the power-law fit and the internal scatter eo is set such that 
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the reduced^ 2 , i.e. the^ 2 per degree of freedom, is unity after 
minimization. <x is chosen to be 200 km s _1 (as in Equation[2]). 
For Equation [3] and H] ctq becomes Mo = 10 n M o and Lq = 
10 n M o , respectively . We use the IDL routine PMFITEXY0 
(Will iams et al.l201Ql) and extend it in order to account for asym- 
metric errorbars. The final values of the fit are given in Table [3] 
which compares the fit parameters from different methods. 



3.2. Treating upper limits 

The methods dev eloped by iFeigelson & Nelson! (Q985) and 
Isobe et alJ (Il986l) combine survival analysis with astronomical 
applications such as determining correlation coefficients and lin- 
ear regression parameters. The basic ingredients are the survival 
function S (z,) (the probability that the object is not detected until 
Zi) and the hazard function A{z/) (the instantaneous rate of detec- 
tion at Zi given that the object is undetected before zi). These 
functions are used to construct statistical quantities and consid- 
ering the non-detections. 

A package of Fortran routines called ASURV0 is available 
for this purpose. The package contains routines for the deter- 
mination of non-parametric correlation coefficients s uch as the 
Cox Proportional Hazard model (Equation 21 in llsobe et all 
1986), the generalized Kendall rank correlation (Equation 28) 
and linear regression methods such as the EM (Expectation- 
Maximization) algorithm with different distributions. We apply 
the routines to our data and list the results in Table |2]for the cor- 
relations and Table |3] for the linear regression parameters. For 
the generalized Kendall rank correlation the routine does not cal- 
culate r itself but the number of standard deviations from zero 
(z - r/[Var(r)]2, where Var(r) is the variance of r under the 
null hypothesis). From this value, the significance of a correla- 
tion between the two variables (P) can be found from a table 
of the integrated Gaussian distribution. In the case of the Cox 
Proportional Hazard model (CPH) a x 1 value is defined in order 
to test for the significance of a correlation and the probability is 
taken from^ 2 distribution tables. 

We stress that survival analysis a) assumes detections to be 
exact measurements, i.e. does not treat uncertainties and b) treats 
upper limits as absolute values, above (or below) where the data 
is definitely not detected. Since neither is the case for our data, 
the results have to be treated with care. 

3.3. Combining uncertainties and upper limits 

In order to treat upper limits and uncertainties in a dataset simul- 
taneously, Kellv (2007) developed a Bayesian method that ac- 
counts for measurement errors in linear regression. This method 
is based on a maximum likelihood approach which is general- 
ized in order to deal with multiple independent variables, non- 
detections (upper limits) and selection effects. 

The basic approach of the routine (available as IDL routine 
LINMIXJ3RR in the Astrolib package) is to generate a random 
sample of the regression parameters drawn from the probabil- 
ity distribution of the parameters given the measured data using 
a Markov chain Monte Carlo algorithm. For each iteration the 
regression parameters and the parameters of the prior density 
are updated until the Markov chain converges to the posterior 
distribution. The saved parameter values can then be treated as 
a random draw from the posterior. The final values are derived 
through the mean or median of the random draw. 



1 Available at: http://purl.org/mike/mpfitexy 

2 Available at: http://astrostatistics.psu.edu/statcodes/asurv 



For performing the Markov chain, two methods are avail- 
able: The first one is the Gibbs sampler, which simulates new 
values of the model parameters and missing data at each iteration 
conditional on the values o f the observed data, current model pa- 
rameters, and missing data Kelly (2007J) . The second method is 
the Metro polis-Hastings algorithm (e.g. [Metropolis et al.l l 1953; 
Hastings 1970) which is used when the selection function is not 
independent of the dependent variables y, the measurement er- 
rors are large or the sample size is small. Since this is the case 
for our sample, we use the Metropolis-Hastings algorithm to de- 
termine the regression parameters (for more details we refer to 
Kellv 2007). In Table[3]we list the parameters obtained with this 
method. For further analysis we adopt the results of the Bayesian 
method as the derived fit. 

4. Discussion 

In this section we discuss the results listed in Table [3] and com- 
pare the resulting correlations of IMBHs with the existing scal- 
ing relations of SMBHs in galaxies. 

4.1. The three major scaling relations 

The first correlations we study are the ones observed in galax- 
ies. Comparing these scaling relations with the ones obtained 
for globular clusters might shed light onto their origin and the 
evolution of globular clusters. Therefore, we first examine the 
relations M. - cr, M. - M, ot , and M. - L. We test for the sig- 
nificance of a correlation by applying the methods described in 
the previous section. Table|2]lists the results of the Kendall t test 
from standard statistics compared to the generalized Kendall r 
and the results from the CPH model from survival analysis. The 
value P gives the probability of an existing correlation between 
the two values and shows high significance (> 90%) for all three 
correlations with all methods. The highest probability is found 
for the M. - cr relation when using the generalized Kendall r, 
but when regarding only the standard statistics, the M. - M tot is 
given the largest significance. It is also the correlation with the 
most stable results over all three methods. 

The existence of these correlations agrees with previous 
measurements of SMBHs in galaxies, where the same major 
scaling relations where found. We apply a linear regression to all 
three relations, using the methods described in the previous sec- 
tions. In Figure|2]we show the three scaling relations, overplotted 
by the best-fits from the Bayesian method (green line in Figure 
[TJ. For the M. - cr relation, all parameters of the linear regres- 
sion from the different methods agree within their errorbars. The 
parameters derived from the different survival analysis methods 
are almost identical. In general, the plots show that the values of 
the survival analysis do not differ much from each other for all 
correlations. The same is true for the FITEXY method and the 
Bayesian approach, which shows that the errorbars have a larger 
effect on the fit than the upper limits. 

Comparing these values with slopes obtained for supermas- 
sive black holes in galaxies (Figure|2] blue line) yields a signifi- 
cant difference for the slope of the correlation between galaxies 
and globular cluster s in the M. - cr re l ation. The most recent 
slope obtained from iMcConnell & M^ (l2012h is (3 a - 5.64 + 
0.32, a factor of two higher than the value that we measured for 
the globular clusters. This is also visible in Figure [2] where the 
values for the galaxies and the globular clusters are overplotted. 
The M. - cr relation seems to be curved upwards when reaching 
the very high-mass black holes and becomes more shallow when 
reaching the low-mass range. 
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Fig. 2. M, - cr, M, - M, ot , and M. - L, , relations of IMBHs and SMBHs in comparison. The slope of the best fit to the GCs (green 
line) is a factor of two smaller than the slope for the SMBHs in galaxies (blue line) for M. - cr, but very similar for the remaining 
two correlations. 



Also for the M. - M tot and M. - L, , correlation all IMBH 
masses lie above the correlation found for galaxies but, as shown 
in Figure [2] the slopes of the correlations do not differ as much 
as they do for the M. - cr relation. In fact, with/? = 1.11 +0.13 
fo r M.-LtM and/3 - 1.05±0 . 1 1 for M. -L tot , the values obtained 
by M cConnell & Ma (2012) agree with the values from Table [3] 
within their uncertainties. 

One possible explanation for the difference in the scaling re- 
lations could be a strong mass-loss of globular clusters in the 
early stages of their evolution. Especially for a> Centauri, NGC 
1851 and Gl this could have had a large effect as thes e objects 
were suggested to be cores of stripped dwarf galaxies (JFr eem an 



ll993MMevlanetalJ200U[Talalietal.ll2rJTllSolhmaet^ 
A higher mass in the early stage of evolution of these objects 
would have shifted the points in all three plots to the right and 
therefore closer to the correlations observed for galaxies. The 
strongest difference in slopes is shown in the M. - cr relation. 
Expansion due to mass loss and relaxation would lead to an de- 
crease of the velocity dispersion and could explain the different 
effect on M, ol and cr. 

Another explanation would be a different mass-radius rela- 
tion for galaxies and globular clusters which complicates the 
comparison of the scaling relations as the velocity dispersion is 
measured from an effective radius in both systems. Since globu- 
lar clusters and galaxies are different objects, formed in different 
environments and p rocesses , this e xplanation seems reasonable 
and is discussed in Gra ham! (1201 lh and references therein. The 
fact that our r elations match with the results found for nuclear 
star clusters ( Graharr feO 1 2t INeumaver & Walchen l2012l) . which 
are found to exhibit shallow slopes, similar to globular clusters, 
supports this theory. Furthermore, clusters in our sample with 
upper limits still have a probability of hosting smaller IMBHs 
not detected up to date. This would bias all the relations towards 
high IMBH masses and could partly explain the difference in 
slopes. In fact, this is indicated by the steeper slopes from the 
survival analysis (where we consider only upper limits but no 
errorbars) as shown in Figure [1] 

4.2. Are there further correlations? 

The three major scaling relations discussed above are known 
for galaxies that host SMBHs at their centers. Even for the 
cases where the slopes of the IMBH correlations agree with the 



SMBH, it does not imply a causal connection. Formation sce- 
narios and co-evolution of host system and black hole are as- 
sumed to be quite different for galaxies and globular clusters. 
For this reason it is crucial to search for fundamental correla- 
tions of cluster-specific values with the mass of the IMBH. This 
has never been done before. This section is dedicated to describe 
the choice of parameters and the results on their correlation sig- 
nificance. 



4.2.1. Distances 

The first parameter we check is the distance to the sun (Dsun)- 
No correlation is expected with this parameter, but it can be used 
to exclude observational biases that might be introduced to the 
measurements. The further away the object, the larger one would 
expect the uncertainty on the black-hole mass to be. Panel a) of 
Figure[3]shows the black-hole mass as a function of distance for 
all Galactic globular clusters in our sample. We exclude Gl in 
the plot (but not in the analysis) as it is the only cluster outside 
the Milky Way and its large distance would distort the figure. 
From the plot no evident correlation between Dsun and M, can 
be found, but the correlation coefficients in Table [2] show a large 
spread reaching from a slight anti-correlation for the standard 
Kendall r to a significant correlation when using survival anal- 
ysis. The reason for this mismatch might arise from the outlier 
position of Gl compared to the other distances. 

As another parameter we study the distance of the cluster 
to the Galactic center (in case of Gl to the center of M31), i.e. 
the position of the cluster in its host galaxy (Dec)- Finding a 
correlation or anticorrelation with galactocentric distance, could 
open speculations about a connection between environment and 
black-hole formation, although the clusters are most probably 
not at the same distance where they were formed. Panel b) and 
the results in Table |2]do not indicate any correlation with Dec- 



4.2.2. Metallicity 

In Panel c) we correlate the black-hole mass with the metallicity 
of the cluster. A cluster with a low metallicity might be more 
likely to form a massive black hole through runaway merging 
of massive stars due to the lack of stro ng stellar winds at low 
metallicities (e.g. IVanbeveren et al.ll2009h . However, according 
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Fig. 3. IMBH masses as a function of several cluster properties such as a) the distance of the cluster to the sun (Dsun), b) the 
distance from the Galactic center (Dec), c ) me metallicity ([Fe/H]), d) the ellipticity (e — 1 - b/a), e) the concentration parameter 
(c), and f) the half-mass radius (ri,) 



to Figure[3]c and the values in Table|2] this does not seem to be 
the case. 



4.2.3. Morphology 

The lower Panels in Figure [3] depict the black-hole mass as a 
function of the three structural parameters: ellipticity (e), con- 
centration (c) and half-mass radius (>/,). The reason why we test 
for correlations in these parameters lies in the possible connec- 
tion between morphology and the cluster being the nucleus of 
an accreted dwarf galaxy. For example, a high ellipticity would 
speak against globular clusters which are assumed to be spher- 
ical and dynamically relaxed systems and for a remnant of a 
dwarf galaxy where aspherical shapes are more common. In the 
same context, an IMBH in a dwarf galaxy is thought to be more 
massive than in a globular cluster since more mass was available 
at its formation. 

Furthermore it has been shown that a black hole prevents 
core collapse and extends the inner regions of the cluster (e.g. 
Baumgardtet al. 2005). Therefore, a correlation of the central 
concentration with black-hole mass would be expected. Also 
clusters with large half-mass radii have been suggested of being 
good candidates of hosting IMBHs at their centers. The highest 
correlation significance is found for the half-mass radius but also 
for the ellipticity a 1 -2 <x detection of correlation is present. For 
the concentration, the values of the different methods show large 
variations ranging from 18% to 67%. In summary, the structural 



Table 3. Best-fit parameters for the three major scaling relations 
obtained form different fitting routines. 





M.-cr 


M. - M T0T 


M. - L T0T 






— FITEXY— 




a 


6.42 ± 0.77 


7.38 ±1.40 


7.96 ± 1.81 


P 


2.21 ± 0.69 


0.69 ± 0.28 


0.76 ± 0.34 


£ 


0.25 


0.33 


0.35 



a 



—SURVIVAL ANALYSIS- 
EM Algorithm 

7.41 ±0.83 8.63 ±2.03 9.84 ±2.10 

3.28 ±0.72 1.01 ±0.34 1.18 ±0.37 

0.36 0.49 0.47 
Buckley- Jones Method 

7.45 8.42 9.47 

3.32 ±0.92 0.96 ±0.32 1.10 ±0.40 

0.42 0.43 0.46 

— BAYESIAN ANALYSIS— 

6.51 ±1.94 7.45 ±3.53 7.49 ± 6.08 
2.34 ±1.63 0.71 ±0.59 0.68 ±1.05 



parameters and the black-hole mass show signs of correlations. 
However, the detection significance is low and needs to be con- 
firmed with a larger data set. 
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5. Summary and Conclusions 

We collected data from the literature from globular clusters that 
were examined for the existence of a possible IMBH at their 
center using stellar dynamics. Our sample consists of 14 globu- 
lar clusters with kinematically measured black-hole masses: six 
of which are thought to host IMBHs, the remaining ones having 
only upper limits on black-hole masses. In order to take uncer- 
tainties and upper limits into account, we use different methods 
to derive the correlation coefficients and linear regression pa- 
rameters. We plot the masses of the central IMBH versus several 
properties of the host clusters, such as total mass, total lumi- 
nosity, and velocity dispersion in order to verify existing corre- 
lations which are found in galaxies, but also with the attempt 
to find new correlations between black-hole mass and cluster 
properties. We find that the three main correlations which are 
observed in galaxies (M, - cr, M. - M tol , and M. - L tot ) also 
hold for IMBHs. The slope of the M. - cr relations differs by a 
factor of two from the fits made to the galaxy sample while the 
remaining correlations are more similar to the ones observed for 
galaxies. 

Furthermore, we test for possible correlations of the IMBH 
mass and other properties of the globular cluster such as dis- 
tance, metallicity, structural properties and size. We find no ev- 
idence for a correlation as strong as the M. - cr relation, but 
we find a trend of black-hole mass with the cluster size, i.e. the 
half-mass radius r/, and its ellipticity. This is reasonable since 
the central black hole leads the cluster to expand by accelerating 
stars in its vicinity. 

We assume in the following that the M. - cr relation is tight, 
has a physical origin and extends to the lowest masses. These as- 
sumptions have not rigorously been demonstrated to date. We as- 
sume further that the IMBHs in globular clusters formed at high 
redshifts and did not evolve much since then, unlike the globular 
clusters themselves which will have suffered mass loss during 
their evolution. With these assumptions, a proposed explanation 
for this behavior is the process of stripping which might have oc- 
curred to these stellar systems when accreted to the Milky Way. 
If they were more massive and luminous in the past they would 
have fit to the correlation of the galaxies instead of being shifted 
to the low-mass end. Especially, these stripping scenarios have 
already been suggested for clusters like Gl and a> Centauri, as 
mentioned in Section [T31 In particular, given the fact that clus- 
ters most probably lost most of their mass during their lifetime 
and galaxies gained mass during their merging history, it would 
be very surprising to find both objects on the same scaling rela- 
tion. 

The fact that the M, - cr shows larger discrepancies when 
comparing to the relation of galaxies than M,-M lol and M.-L tot 
might arise in either different mass-radius relations of these very 
different black-hole host systems as well as the expansion driven 
decrease of the velocity dispersion in globular clusters. This is 
supported by the fact that the scaling relations agree with the 
upper limits for nuclear clusters, stellar systems more similar to 
globular clusters. N-body simulations and semi-analytic expan- 
sion models are needed in order to quantitatively study these the- 
ories. We note that also a bias towards higher black-hole masses 
due to detection limits can not be excluded. 

For the future it is desired to extend this sample and con- 
firm the reported correlations and their slopes. With integral- 
field units combined with adaptive optics, the search for IMBHs 
can be extended to extragalactic sources. Nevertheless, a com- 
plete search of Galactic globular clusters will provide a larger 
sample and help to further constrain critical observables which 



hint towards the presence of an IMBH in the center of a globular 
cluster. 
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